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A  MODEL  IN  WHICH  COMPONENT  FAILURE  RATES  DEPEND  ON  THE  WORKING  SET 


by 

Sheldon  M.  Ross 


1.  INTRODUCTION 

Consider  an  n  component  system  having  some  arbitrary  monotone  coherent 
structure  (see  Barlow  and  Proschan  [1]  for  suitable  definitions).  We  sup¬ 
pose  that  each  component  Is  Initially  on  and  stays  on  for  a  random  time  at 
which  It  fails.  The  problem  of  Interest  is  to  characterize  the  distribution 
of  the  time  until  the  system  falls.  Whereas  this  problem  is  usually  con¬ 
sidered  under  the  assumption  that  the  component  lives  are  independent,  we 
are  concerned  with  the  following  model  which  allows  for  dependencies  in 
these  life  distributions:  We  suppose  a  Markovian  model  in  which  the  failure 
rate  of  a  given  component  at  any  time  is  allowed  to  depend  on  the  set  of 
working  components  at  that  time.  Specif ically, we  suppose  that  if  at  some 
time  W  ,  W  C  {1,2,  ...,  n}  ,  represents  the  set  of  working  components  then 
for  i  e  W  the  instantaneous  failure  rate  for  component  i  is  X^(W)  . 

We  start  by  giving  a  sufficient  condition  for  the  distribution  of  system 
life  to  be  NBU  where  we  say  that  the  nonnegative  random  variable  T  has  a 
NBU  (new  better  than  used)  distribution  if 

P{T  >  a  +  t  |  T  >  s}  <  P{T  >  t)  for  all  s,t  >.  0  . 

In  words,  the  above  states  that  the  probability  a  used  item  survives  an 
additional  t  time  units  is  less  than  the  corresponding  probability  of  a 
new  item.  We  are  now  ready  to  show  that  if  the  failure  rate  of  a  component 
increases  as  the  set  of  working  components  decrease  then  system  life  is 


NBU. 
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Proposition  1; 

If  for  all  sets  , 

X^Wj)  >  X±(W2)  ,  i  e  Wx  , 


then  the  time  until  system  failure  is  NBU. 

Proof: 

For  any  set  of  components  U  ,  let  T^  denote  the  time  until  the 
system  fails  when  W  consists  of  the  set  of  components  that  are  initially 
working.  He  will  start  by  showing  that  if  Z  C  W  then  T_  <  T  .  The 

u  "T  W 
St 

proof  of  this  will  be  by  Induction  on  k  -  |z|  +  |w|  ,  where  |u|  equals 
the  number  of  elements  in  U  .  It  is  obvious  for  k  -  1  (for  in  this 
case  Z  «  4>  and  so  Tz  =  0) ,  and  so  assume  it  whenever  |z|  +  |w|  -  k  . 
Now  suppose  that  Z  CM  and  |z|  +  |tf j  -  k  +  1  .  For  i  e  Z  ,  define 
X^  to  be  art  exponential  random  variable  with  rate  X^(Z)  .  Also  for 
J  e  W  -  Z  define  Yj  to  be  exponential  with  rate  X^  (W)  .  In  addition, 
suppose  that  all  the  X^  and  Y^  so  defined  are  Independent.  Now  let 

X  -  min  | min  X.  ,  min  Y  . 

(ieZ  jeW-Z 

There  are  two  cases  we  need  consider: 

Case  1:  X  ■  X^  for  i  e  Z 


In  this  case,  we  can  set 
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X  +  T 


{Z-i} 


*  Xi(W) 

lX  +  T*W-i}  with  probability  ^ 


TW"' 


X  +  T, 


Ai(W) 


W 


with  probability  1  -  ~(z) 


where  T*  is  meant  to  be  a  random  variable  Independent  of  all  the 

and  Yj  and  with  the  same  distribution  of  Ty  .  Since  {Z  -  1}  C  {W  -  i}  C 

H  ,  it  follows  by  the  induction  hypothesis  that  in  this  case  T„  <_  TL,  . 

Z  st  W 


Case  2:  X  »  for  some  j  e  W  -  Z 
In  this  ca8e,we  set 


tz  -  x  + 

’v  "  x  +  T{w-j}  • 

j}  ,  it  again  follows  by  the  Induction  hypothesis  that 

Hence,  for  Z  C  W  ,  T_  <  Tu  .  Now  suppose  all  components  are  Initially 

1  st  W 

on  and  that  the  system  is  still  working  at  time  s  .  Now  no  matter  what  the 
set  of  working  components  is  at  time  s  ,  it  follows  from  the  above  that  the 
remaining  life  is  stochastically  smaller  that  T,-  «  ,  which  proves 

the  proposition. | | 


As  Z  C  {W  - 
Z 


T-  i  Tw 


st 


Regards: 


(1)  Proposition  1  need  not  be  true  without  the  monotonicity  assumption 
on  A±(W)  .  For  a  counterexample,  consider  a  parallel  system  with 
n  *  2  and 


X1(l,2)  -  A2(l,2)  -  1 
X1(l)  -  e  ,  X2(l)  -  1  . 

Suppose  the  system  is  working  at  t  .  Now  as  t  becomes  larger  at 
some  point  the  system's  failure  rate  starts  to  decrease  because  it 
becomes  more  and  more  likely  that  only  component  1  is  working  (the 
only  other  possibility  of  any  probability  being  that  only  2  is 
working).  Hence,  system  life  will  not  be  NBU. 

(ii)  As  the  failure  rate  of  a  working  component  depends  on  the  set  of 

failed  components,  the  question  arises  as  to  whether  Proposition  1 
would  remain  true  if  this  failure  rate  were  allowed  to  depend  on 
the  order  in  which  the  components  have  failed.  That  is,  suppose 
that  X^dj^ij,  ...»  i^)  ,  i  ^  ij  ,  J  ■  1,  ...,  k  ,  is  the  failure 
rate  of  component  i  when  components  i^,  ...,  i^  have  failed  and 
in  that  order.  Would  system  life  be  NBU  if  X^(i^,  ...,  i^)  <_ 
1^(1^,  ...,  i^.ij^)  ?  The  answer  is  no  for  consider  the  following 
example  for  a  parallel  system: 

n  ■  3  ,  X:  -  X2  -  10^  ,  X3  «  1 


A. (1,2)  -  10 


X.(2,l)  -  1 
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where  X  is  the  initial  failure  rate  of  component  i  .  How  after 
a  short  time  both  1  and  2  would  be  failed  and  so  the  remaining  life 
will  be  a  mixture  of  an  exponential  with  rate  10  (if  1  failed  before 
2)  or  an  exponential  with  rate  1  (if  2  failed  before  1) .  But  a 
mixture  of  exponentials  with  unequal  rates  has  a  decreasing  failure 
rate  and  so  system  life  could  not  be  NBU. 

(iii)  When  n  »  2  ,  the  joint  distribution  of  the  lifetimes  of  the  two 

components  is  called  the  Freund  Distribution  (see  [2]).  It  can  be 
shown  in  this  case  (see  [6])  that  the  time  of  system  failure  has  the 
stronger  than  NBU  property  of  being  an  increasing  failure  rate  on 
average  (IFRA)  distribution.  We  do  not  know  if  this  result  can  be 
extended  to  the  case  n  >  2  . 


(lv)  An  interesting  special  case  obtains  when  we  take 


X±(W) 


xlC 


1 

l  ‘j 

jeW  J 


i  e  W 


where  a^  are  given  nonnegative  numbers.  Such  a  situation  would 

arise  from  the  following  weighted  load  sharing  model:  Suppose  that 

an  n  component  system  is  subject  to  a  constant  load  pressure  C 

which  must  be  allocated  among  the  working  components .  Suppose  also 

that  the  allocation  is  determined  by  a  set  of  weights  a^,  ...,  aQ 

such  that  if  at  any  time  W  is  the  set  of  working  components  then 

the  load  taken  on  by  component  i  ,  i  e  W  ,  is  C  a. /  1  a.  .  If 

1  jeW  J 

in  addition  we  suppose  that  the  failure  rate  of  component  i  is 
proportional  to  (with  proportionality  constant  X^>  the  load,  it  is 
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as aiming  then  the  ebove  obtains.  For  this  model, it  can  be  shown 
(see  [5])  that  for  a  parallel  system  the  time  until  system  failure 
(that  is  the  time  until  all  components  have  failed)  is  an  increasing 
failure  rate  (IFR)  random  variable. 


« 


( 


*• 


2.  LAPLACE  TRANSFORM  OF  SYSTEM  LIFE 


For  a  given  structure  one  can,  upon  conditioning  on  the  order  in  which 
the  components  fail,  obtain  an  expression  for  the  Laplace  transform  of  T  , 
the  time  of  system  failure.  For  instance,  suppose  a  parallel  system  which 
fails  when  all  components  fail.  Then  letting 

A^(W)  *  Xi(WC)  ,  where  Wc  ■  complement  of  W 


we  have 


E[e'sT]  - 


l 

(i^». . .  * in) eJP 


F(i^,  • • • i  iQ) 


n 

n  X  *i,(il*  ik-l) 

n  - 

»+  I  X.  (i.»  ...»  i.  ,) 

L  J-k  J 


where 


£  is  the  set  of  all  nl  permutations  of  1,2,  ...,  n 


and 


^*1'*2'  *  *  *  *  "  Q 


Xi2(1l) 


...  1 


I  X.  I  X.  (i.)  I  X,  (i^*  *•»*  ij^j) 
j-1  *J  J-2  li  1  J-k  1  ^ 


is  the  probability  that  components  fall  in  that  order. 

The  above  can  easily  be  understood  by  noting  that  given  that  the 
components  fail  in  order  ...»  iQ  the  time  between  successive 
failure  components  are  independent  exponentials  with  rates 
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3.  SIMULATING  THE  PROCESS 

Let  denote  the  failure  time  of  component  i  .  The  random  vector 

(Tp  ...»  Tn>  can  most  easily  be  simulated  as  follows:  Let  , 

i  ■  1 . n  be  independent  exponentials  with  respective  rates  , 

i  *  1*  ...»  n  .  Now  order  the  X.  and  let  i  be  such  that 

i  1  /  n 


n 


The  above  follows  by  first  noting  that  given  ip  iQ  ,  it  follows  from 

the  lack  of  memory  property  of  the  exponential  that 


are  independent  exponentials  with  rates  1.  The  denominator  term 
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4-t,  *ij(il*  '*•»  ik-l)  in  the  definition  of  T  thus  gives  the  exponential 

J-k 

its  appropriate  rate. 


Remark: 

When  n  -  2  and  ^  ^(j)  i  j  M  ,  the  above  expresses  T2,T2  38  an 
increasing  homogeneous  function  of  X1,X2  .  This  was  noted  and  used  in  [6] 
to  show  that  the  are  associated  (follows  from  the  fact  they  are  in¬ 

creasing  functions  of  the  independent  random  variables  X±  (see  [1]  for  a 
proof  of  this))  and  also  that  system  life  is  IFRA  (follows  from  the  fact  that 
the  function  is  not  only  increasing  but  also  homogeneous — see  [4]  for  a  proof 
of  this).  Unfortunately,  when  n  >  2  ,  these  functions  are  no  longer  in¬ 
creasing  . 


10 


4.  THE  MODEL  WITH  REPAIR 

Let  us  suppose  that  failed  components  are  repaired.  Specifically, 
suppose  that  the  repair  rate  of  component  i  when  then  set  of  working 
components  is  W  ,  i  i  W  ,  is  u^(W)  .  This  gives  rise  to  a  continuous 
time  Markov  chain  with  2n  states — all  possible  subsets  W  C  {1,2,  ...,  n}  . 
To  solve  for  the  steady  state  probabilities  is  in  general  a  difficult  task 
but  it  simplifies  in  the  following  special  case: 

Special  case: 

For  functions  f(k)  ,  g(k)  ,  k  ■  0,1,  ...,  n  and  positive  constants 
^k  *  **k  ’  ^  •  •  • ,  n  , 

XjCw)  -  X  f(|w|)  ,  j  e  W 
M±(W)  -  U^dwj)  ,  i  i  W  . 


Proposition  2: 

Under  the  conditions  of  the  above  special  case,  the  stationary  prob¬ 
abilities  of  the  set  of  working  components  is  given  as  follows:  For 
W  -  Uj^,  ....  ij^J  » 


Wi  ^i  •  •  * 

P({i  4  })  ,  1  2  \  g(k-  1)  ...  g(0) 

•••’  V'  \4  x4  ...  x.  f(k)  ...  f(i) 


P(*> 


ll  2 


where  P($)  is  the  stationary  probability  that  all  components  are  failed 
and  can  be  obtained  by  summing  the  above  over  all  W  and  equating  to  1. 
In  addition  the  chain,  in  steady  state,  is  time  reversible. 
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Proof : 

To  verify  the  above,  all  we  need  check  Is  that  the  proposed  stationary 
probabilities  satisfy  the  time  reversibility  equations.  That  Is  we  need 
check  that,  for  the  proposed  stationary  probabilities, 

PU^ . ijJU^ftt)  -  P«i2 . i^})^  8(k  ~  D  • 

But  this  is  iamediate  and  so  the  result  follows. | | 
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